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Abstract 

Finite length of a one channel wire results in crossover from a Tomonaga- 
Luttinger to Fermi liquid behavior with lowering energy scale. In condition 
that voltage drop (V) mostly occurs across a tunnel barrier inside the wire 
we found coefficients of temperature/voltage expansion of low energy conduc- 
tance as a function of constant of interaction, right and left traversal times. 
At higher voltage the finite length contribution exhibits oscillations related to 
both traversal times and becomes a slowly decaying correction to the scale- 
invariant y 1 / 3-1 dependence of the conductance. 
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Quantum transport in Tomonaga-Luttinger liquids (TLL) has attracted a great deal of 
interest as it was suggested to be realized in a ID constriction and the edge state of the 
fractional quantum Hall (FQH) liquid |3|]. Both suggestions have been supported in recent 
experiments on FQH and on ID H transport. Even more experimentalists |J made claim 
on observation of the interaction effects in the ID transport, however, without comprehen- 
sible connection with a theoretical model. Unless there is a resonant tunneling the repulsive 
interaction typically suppresses conductance at low energy. The experiment by Tarucha's 
group H on ID transport through a long wire with a weak impurity random potential inside 
demonstrated a crossover from TL liquid to Fermi liquid behavior at low temperature. This 
crossover is a finite length effect |7j and may be described in an inhomogeneous TLL model 
(ITLL) |8|-[10|. The ITLL model predicts both the conductance behavior in the Fermi liquid 



region up to a renormalization constant (i. e., relations between the coefficients at tem- 
perature/voltage in different integer degrees) and the interaction dependent non-analytical 
behavior in the TLL region. 

The aim of this work is to examine this crossover in the low voltage conductance of the 
one channel wire where its suppression is mostly determined by a high point barrier located 
inside the wire. Position and height of the barrier are assumed to be due to an external 
gate. Therefore both distances Lr(l) from the barrier to the right/left reservoir and the 
ratio ( = Lr/(Lr + L&) are assumed to be known. Then it will be shown below that 
just two more parameters, which are the total traversal time equal to sum of the right 
and left ones:t = t R + t L and the constant of the forward scattering g, are necessary for 
desciption of the crossover, and that the ITTL model gives quite a few ways to determine 
these parameters from either high voltage (V > 1/to) or low voltage (V < 1/to) conductance 
measurements. Temperature dependence of the conductance in this model was considered 



in 



TTJ. Similar model turned out to be useful for experimental study of the FQH transport 
It was noticed recently |12] that there is a special way of connection between a v — 1/3 
FQH liquid and leads when the whole setup corresponds to the g = 1/3 ITLL model. If so, 
the results obtained below for spinless electrons could be directly addressed to that FQH 



device when the FQH liquid embeds a point scatterer. 

Under condition that the link between two parts of the wire is weak enough it suffices 
to apply the tunneling Hamiltonian approach in the lowest order to describe the transport 
P JIl||l4jj . Current flowing through the weak link located, say, at x = inside the wire is 
given by the operator J (t) = —i[Aif)R(0,t)if)L(0,t) — h.c], (e,h = 1), where A is the tunneling 
amplitude and iI>r,l{x, t) are the electron annihilation operators in the right < x < Lr and 
in the left — Ll < x < part of the wire, respectively. The average current under voltage V 
applied to the left lead can be written as: < J >= 2tt\A\ 2 J de[f(e—V) — f(e)]pR(e)pL(e—V), 
where / is Fermi distribution. The problem reduces to finding of the tunneling density of 
states of the right (left) end of the junction Pr(l)(c) which are the sum of the particle 
p p {e) = (1 — /(e))p(e) and hole p^(e) = /(e)p(e) densities. The latters relate to the particle 
correlator as p p (e) = l/(27r) / dte let < i[)(0, t)ip + (0, 0) > and to the hole one as Pft(e) = 
1/(2tt) J dte iet < ^+(0, 0)^(0, t) >. 

Tunneling density of states - To calculate the tunneling density of states pR on the right 
side of the weak link let us first consider spinless fermions and apply bosonization to the ip 
field under condition of an elastic reflection from the boundary located at x = [r5,IB|. ( 



Carrying out this calculations we will omit index "R" below. ) Bosonic repersentation of 
the ip field reads ip(x,t) = J2 a =r,i 4>a(x,t) = (2iia)~ 1 J2± ex P{i(@( x ,t) ± <!>{%, t))/2}, where 
ip r {i) is the right (left) going chiral component of ip and the 9 and 4> fields are bosonic and 
mutually conjugated [9(x,t),(j)(y,t)] = 2nisgn(x — y). The elastic reflection means that 
ipi(0,t) = e lS ipr(0,t) with an appropiate phase shift 5. This results in both: 

0(0, t) = 5, ^-d x 9(x, t) \ x=0 = </v+(0, t)W(0, t) - ^ + (0, t)V/(0, t) = 0. (1) 

Then the density of particle states could be found as 

p U) = [ +0 ° dte iet+±[<e(o,t)e m >-<e 2 m>] (2) 

2n J-oo 

where the value of the free electron tuneling density was introduced as: po = (1+cos 5) / (irv). 
The problem reduces to finding the 9 field correlator. It can be done for the finite length 



piece of the wire adiabatically connected to the lead making use of the ITTL model [pHlCfl . In 



this model the Tomonaga-Luttinger interaction {J2 r ,i Pa) 2 is switched on in the wire x < La 
and switched off outside. Then the Hamiltonian takes a bosonized form 



where function ensuing from the interaction can be approximated in the low energy 

limit by a step-function: u(x) — 1 if x > Lr and u(x) — u — 1/g < 1, otherwise. The 

correlator of the 9 field ordered in imaginary time T(x,y,r) =< T T 6(x,T)6(y,0) > can be 

shown to satisfy the following equation 

1 Att 

{-TYf-A + 9l}T(x, y, r) = S(x - y)S(r) (A) 

v z u z (x) V 

under the boundary conditions d x T(x, y, t)\ x= q = following from (|l|). Fourier transform 

of this correlator T(x,y,u) is symmetrical under uj — > —uo. It can be constructed from the 

solutions of the homogeneous equation corresponding to Eq.(^) 

{ -TTTs ~ d DU x ) = i-T^ ~ ^}M*) = (5) 
v 2 u 2 (x) v 2 u 2 (x) 

At: 

T(x,y,u) = [0(x - y)fu{x)K(y) + 6(y - x) f w {y)K{x)} 

if these solutions meet boundary conditions: /^(0) = 0, fw(x) = exp(—uix/v) at x — ► 00 
and positive ui. The Wronskian W{ui) is equal to — /^(0)/i^(0) and, hence, T(0, 0, u) = 
— Att/[v (lnf u (0)Y]. The only solution we need can be written as right going plus reflected 
left going waves at x < Lr. The reflection amplitude r = r# = —e~ 2ri , (tanh(ry) = 1/u) 
for the 9 field is negative for the repulsive interaction. It is related to the one for the 
field |5|,[nj as r 9 = ~ r 4> by the duality symmetry. Substituting this solution one can find 
T(0, 0, to) = — t&nh(utji + r}) with tji equal to the time of travelling from the junction to the 
right lead. Analytical continuation of this function [— T(0, 0, —iu) + 0)} is the retarded Green 
function for the 6 field. Imaginary part of the latter multiplied by the Bose distribution 
function for holes 1 + /s(tf ) and by a factor (—2) coincides with the Fourier transform of the 
correlator at to. Then the particle density of states @ in dimensionless units is obtained as 

P P (e) = / dpexp{iep + / du^H 1 + / fl M ~ , b S > ( 6 ) 

/7T7 J-00 Joo tanh(77j 



where the inverse temperature ft and energy e,co were scaled as ft = l/(Ttu), e = t^e, uj = 
tnuj (T is the temperature in energy units). The dimensionless cut-off parameter 7 is 
(EptR,) -1 , and p is dimensionless time p = t/tR. The hole density of states Ph{s) can 
be found as Ph(s) = Pp(—e)- 

The correlator used in ([]) could be represented as a product after expansion of the 
imaginary part of the tangent in a sum of exponents: 

E F ( -ij-K/ft \ u ^ ( |sinh((2n + i 7 )7r//3)| 2 \ ur " 



< W(0,J#7(0,0) >= 



2irv Vsinh((p - i^/ft) JJ \j[ ± s inh((2n ± (p - i-y))ir/ft) t 



(7) 



This expression can be easily understood as the product of the contributions of the 2n 
length paths connecting (0,p) and (0, 0) points and undergoing n reflections from a x = Lr 
non-elastic boundary with the negative reflection amplitude r = — e~ 2v and n reflections 
from the x = elastic boundary with unit reflection amplitude. The product has a good 
convergence due to exponential decrease of r n with n. Substituting it into one can 
find low (e <C 1) and high (e ^> 1) energy behavior of the tunneling density of states. 
Similar calculations in the spinful case require change of Jm[tan(a) + if])}/ tanh(r/) in (^|) 
into (1 + Jm[tan(d) + if])}/ tanh(^))/2. The product representation of the correlator can be 
obtained from (I7|) if the exponent of the second multiplier u is changed into (u + l)/2 and 
the exponents of all the rest of the product ur n into ur n /2. 

The calculation of the low energy expansion is tiresome but straightforward. Up to the 
fourth order in energy it gives in the spinless and spinful cases, respectively: 



Pie) 



Pie) 



c lPo i u - l [l + {l-g 2 " 



0.104e 2 (^j +0.01e 4 



0.104e 2 (Jj +0.01e 4 



2 + 6 [ft 



I - g 2 )(4.5 - Q.5g 2 ) 
+ (1-/)[1.34-1V] 

+ (l-2 2 )(2.22 - 3.22(7 2 ) 

+ (1-g 2 ) [0.644 - 0.926<? 2 ] (j) ) 



(8) 



T + 12 1/3 



(9) 



where the constant c% = n n> o{2n) 2urn modifies a renormalization parameter y _1 and is 
model dependent. Except for this parameter the functions p(e) and, hence, all coefficients 
of (|8]|D in the curly brackets are universal functions of g and £r. The coefficients in the 
spinful case are approximately two times less than the ones in the spinless. In the second 
order this correspondence is exact. 

The high energy behavior of the tunneling density at zero temperature and positive e 
can be obtained making transform of the counter of integration in @ into a sum of the 
contours going around the cuts of the correlator (|7|) in the complex time plane from ioo + 2n 
to 2n and back. The result takes the form: p(e)/po = — — (sin(7ra)r(l — u)^"^ 1 + 2r(e)) 
in the spinless case. The spinful expression just needs replacement of u by (u + l)/2 in the 
spinless one. It reveals that finite length corrections 

r(e) = a n sm{irur n )r{l - ur n ) cos(2ne - ^(1 - r"))^'"" 1 (spinless) (10) 

n>0 ^ 

r(e) = Y, bn sin(7rur"/2)r(l - ur n /2) cos(2ne - -(2 - r «))^"/2-i ( spinM ) (n) 

n>0 ^ 

to the scale-invariant power law dependence display an interference structure. Here 
the constant a n in the spinless density of states flT0P is given by a" 1 = 
2u n u(i-r n ) j] m _^ n>0 l — (iLj and the constant b n in the spinful density (|TTj) is 6 n = 



a n / (2n). Both of them are on the order of 1 at small n and decrease with increase of the 



number. However, a quick convergence in (JTO, 11) is mostly due to the ur n factor. The first 



two terms of the sums bring the leading contribution. They describe weakly decaying oscil- 
lations with the periods equal to 2t R (n = 1) and 4i R (n = 2), respectively. The second term 
always dominates at large energies as r < 0. However, the crossover energy is exponentially 
large at small u. Comparing coefficients of the first two oscillating terms of the sum in ( pT|) 
one can gather that the finite length oscillations of the 2tn period dominate the ones of the 
Atn period above £ ~ 1 unless u > 1.8 in the spinless case or u > 2.1 in the spinful one. It 
is in accordance with the numerical calculations ||17|| . This structure in the density of states 
may be understood as a quantization of the plasmonic modes inside the wire [|T8]]. 

Conductance - It is convenient below to redefine the dimensionless energy in units of the 



inverse total traversal time as T = toT, V = toV. Low energy conductance is an analytical 
function of T, V < (1 — C) _1 - Integrating the low energy tunneling density of states 
(§,0), it can be found in the spinless and spinful cases, respectively: 



a 



c 2^2(„-l) 



R0C1 



u-l 



i +(i-s 2 R 2 



v 2 W 

2 3 



+ 



-Ci + (1-0 , 
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-At<7 



a 



C17*- 1 



(u-l)/2 



1 + 



V 2 (ttT)- 

~2~ + 3 



+ 



2N2 



6 64 



+ (1-5 



2 ,9-13^ 



32 



ttTV) 2 + A T a 



{■kTV) 2 
(12) 



(13) 



where 7 _1 = to-^F and Ci = £(1 — 0? C2 = 1 — 2Ci,C3 = C 4 + (1 — C) 4 are geometrical 
coefficients determined by the only parameter £ of the barrier position and Rq 1 is a free 
electron conductance of the junction. At<j is produced by the finite difference between the 
temperatures of the right and left reservoirs AT = Tr — Ti and T = (Tr + T^)/2. In 
the lowest linear order it is non-zero unless ( = 1/2 and reads for the spinless and spinful 
conductance, respectively: 

c ^2(«-l) (1 _^ )(9 _ 13 ^ 



At<j 
At& 



Rocr 1 72 

i-g 2 



-UWtat 



C17"" 1 



(0.74 - 1.073£ 2 )C 4 WTAT 



(14) 
(15) 



RoC (u~m 12 

where (4 = £ 4 — (1 — Q 4 . For = 1/2 Aj-a in([l4D, (p~5|) vanishes and the conductance is 
always increasing due to (AT) 2 term as the temperature difference is introduced. Relation 
of the T 2 term coefficient to the V 2 one is universal. Following Weiss's consideration [EC 



one 



can prove it is equal to (27r 2 )/3 in all orders of the perturbation expansion in the tunneling 
amplitude [[HJ. This is similar to Wilson relation in the Kondo problem. One can see from 
(|j~2|), (|i~3| ) that both the interaction constant g and the traversal time t may be determined 
from the low energy expansion by comparing coefficients for the terms of the second and 
fourth order in energy. 

Next we consider the high energy behavior of the conductance. To describe the lat- 
ter for V ^> 1/(1 — (),!/( it is convenient to start from a representation of the current 



in the form: < J >= \A\ 2 J dt{e lVt - e~ iVt ) < ^(0,t)^(0,0) >< V£(°> t)i/> L {0, 0) >. 
The correlators may be expressed in the form (|7|). Making asymptotic integration we 
come to the conductance which in the spinless case could be written as dJ/dV = 
(tt RoY 1 (toE F ) 2i ~ l ~^ [sin tt(2u - 1)T(2 - 2u)V 2u - 2 + Aa(V)}. Expression for the spinful con- 
ductance is obtained by replacing u by (u + l)/2. At zero temperature the finite length 
corrections Aa(V) of spinless and spinful conductance may be found for irrational values of 
(, respectively, as 



Aa(V) = 2V 2u ~ 2 



n>0 



cos(2nVC + tyn,R) cos(2nF(l - C) + 2^ L ) 



( y C) 2«(l-r"/2)-l^ (_C_) {V (1 _ C))2u( l- r «/2)-l^ ( i_C) . 

x a n sin(7TMr n )r(l - ur n ) (16) 



Aa(V) = 2V U ~ 1 

n>0 



sin(2nVC + ip n , R ) sin(2nV(l - () + ip niL ) 



x b n sin(vTMr n /2)r(l - ur n /2) (17) 



where g^(x) = {2n) u ~ 1 IT m >o 1 ~~ (xn/m) . The phase shifts are (r" — (u — l)r^ nx ')n / 4 
([nx] denotes the integer part of nx) with x — — C) f° r <Pn,R and x = (1 — £)/C for <^n,L- 
Again one can see that only the n = 2 term is important at high V. However, if — r <C 1 
the n = 1 oscillations dominate over a large range of energy above 1. With decrease of 
C < 1/2 the oscillations of the 7r/(£to) periodicity acqure much larger amplitude than those 

Of 7T/(1-C). 

If ( is rational: C/(l — C) = ^2/^1 the resonant enhancement of the oscillations with the 
frequency 2n\C,V = 2n 2 (l — ()V occurs. Since only lowest n's contribute to the sum (|I~7D 
this resonance is important when both n^'s are small: ( = 1/2,1/3,2/3. For ( = 1/2 the 
finite length correction of the spinless and spinful conductances reads 



Aa(V) = V 2u ~ 2 4na 2 n cos{2nVC + vr(r" - 1)) sm(2nur n )T(l - 2w n )(\/C) Mr "~ 1M (18) 

n>0 

Aa(V) = V 11 - 1 2a n sin(2^C + nr n /2) sin(7rar")r(l - ur^VQ 11 ^-^ (19) 

n>0 

The resonant enhancement strengthens oscillations of the 4( frequency and weakens the ones 
of the 2(. These interference structures are shown in Fig.l for spinful electrons. 



S 



With decrease of V the high voltage behavior of the conductance (|16D, ([17]) will meet 
the low voltage one ([12]), (|13|). However, if ( <C 1/2 there is a region l/(l-()<cV<l/( 
where conductance relates to the low and high energy tunneling densities of states of the 
right and left parts of the wire, appropriately. In the leading order it is described by 
a = p R (0)p((l-()V)/(R o p 2 o). 

In summary, we considered suppression of the conductance through a tunneling barrier 
in a ID channel constriction in the ITLL model accounting for the finite length of the wire. 
Starting from the zero energy the conductance increases analytically as T 2 , V 2 and higher 
even integer degrees. The relations between all coefficients are determined by the constant 
of interaction and two traversal times. At voltages higher than inverse of the sum of the 
traversal times conductance approaches the infinite length scale-invariant dependence. The 
deviation from the latter produced by the finite length contribution decays as a negative 
degree of the voltage oscillating with the periodicities related to double traversal times for 
a weak repulsion and to quadrupole ones if the interaction become stronger. 

We acknowledge S. Tarucha for useful discussions. It is a special pleasure for one of us 
(V.P.) to thank T. Iitaka for his help in conducting calculations. This work was supported 
by the Center of Excellence and partially by the fund of the JSPS for development of 
collaboration between the former Soviet Union and Japan. 
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FIGURES 



Dependences of the finite length correction (EptoY 1 '^ Aa(V) / Ro to the differential con- 
ductance on volage V measured in n over traversal time t unit for spinful electrons at 
u = 1.396: the solid line relates to the symmetrical position of the weak link; the dotted 
line to the case when the ratio between the lengths of the right and left shoulders is 1/4. 



10 



REFERENCES 
[1] W.Apel and T.M.Rice, Phys. Rev. B 26, 7063 (1982). 

[2] C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 1220 (1992); Phys. Rev. B 46, 
15233 (1992). 

[3] X.G.Wen, Phys.Rev. B 40, 7387 (1989); K.Moon et al , Phys. Rev. Lett. 71, 4386 
(1993). 

[4] A.M. Chang, L.N. Pfeiffer and K.W. West Phys. Rev. Lett. 77, 2538 (1996); F.P 
Milliken, C.P.Umbach, and R.A. Webb, Solid State Commun. 97, 309 (1996). 

[5] S. Tarucha, T. Honda, and T. Saku, Solid State Commun. 94, 413 (1995). 

[6] A. Yacobi et al, Solid State Commun. 101, 77 (1997); K.J. Thomas et al, Phys. Rev. 
Lett. 77, 135 (1996). 

[7] M. Ogata and H. Fukuyama, Phys. Rev. Lett. 73, 468 (1994). 

[8] D. L. Maslov, Phys. Rev. B 52, R14368 (1995); D. L. Maslov and M. Stone, Phys. Rev. 
B 52, R5539 (1995). 

[9] V. V. Ponomarenko, Phys. Rev. B 52, R8666 (1995); Phys. Rev. B 54,10328 (1996). 
[10] I. San and H. J. Schulz, Phys. Rev. B 52, R17040 (1995). 
[11] A.Furusaki and N.Nagaosa, Phys. Rev. B 54, R5239 (1996). 
[12] D. Chklovshii and B. Halperin, |cond-mat/9612127 . 
[13] F. Guinea et al , Europhys. Lett. 30, 561 (1995). 

[14] PFendley, A.W.W.Ludwig and H.Saleur, Phys. Rev. Lett. 75, 2196 (1995), Phys. Rev. 
B 52, 8934 (1995). 



11 



[15] J.L.Cardy, Nucl. Phys. B275[FS17], 200 (1986); Nucl. Phys. B324, 581 (1989); 
H.W.J.Blote,J.L.Cardy and M. P. Nightingale, Phys. Rev. Lett. 56, 742 (1986); 
M.Ameduri, R.Konik and A.LeClair, |hep-th/9503088; . 



[16] S.Eggert and I.Affleck, Phys. Rev. B 46, 10866 (1992); M.Fabrizio and A.O.Gogolin, 
Phys. Rev. B 51, 17827 (1995); S. Eggert et al , Phys. Rev. Lett. 76, 1505 (1996). 



[17] V.V. Ponomarenko, |cond-mat / 9703085 



[18] Y.V.Nazarov, A.A.Odintsov and D.V.Averin, Europhys. Lett. 37, 213 (1997). 
[19] V.V. Ponomarenko and N.Nagaosa, in preparation. 
[20] U. Weiss, Solid State Commun. 100, 281 (1996). 



12 



V. Ponomarenko and N. Nagaosa 



Ac 




Fig. 1 



